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2 Final Merger of Black-Hole Binaries 

Abstract Recent breakthroughs in the field of numerical relativity have led to dramatic progress 
in understanding the predictions of General Relativity for the dynamical interactions of two black 
holes in the regime of very strong gravitational fields. Such black-hole binaries are important as- 
trophysical systems and are a key target of current and developing gravitational-wave detectors. 
The waveform signature of strong gravitational radiation emitted as the black holes fall together 
and merge provides a clear observable record of the process. After decades of slow progress, 
these mergers and the gravitational-wave signals they generate can now be routinely calculated 
using the methods of numerical relativity. We review recent advances in understanding the 
predicted physics of events and the consequent radiation, and discuss some of the impacts this 
new knowledge is having in various areas of astrophysics. 
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1 Introduction 

1.0.1 Black Holes It has been nearly a century since Albert Einstein's 
profound physical insight revealed our current standard model of gravitational 
physics, General Relativity. Among the theory's extraordinary consequences was 
the predicted existence of black holes, nonlinearly self-gravitating, stable objects 
in which gravitational forces have completely overcome all other physical inter- 
actions. Though once considered a mathematical curiosity, General relativity's 
description of black holes now provides the best explanation for a widespread 
class of astrophysical objects whose gravitational potential wells power many of 
the most energetic astronomical phenomena observed. These range from stellar 
black-holes powering X-ray sources in the neighboring regions of our galaxy, to 
supermassive monsters with masses ~ (10 6 - 1O 9 )M , where M Q is the mass of 
the Sun, with far-reaching astrophysical consequences. 

1.0.2 Gravitational Wave Observations In the coming decade, antic- 
ipated observations of gravitational waves from black hole binaries will open a 
new window onto the universe. Interpreting such observations will richly engage 
aspects of our theoretical understanding of strong-field gravity which have never 
before been confronted with empirical observations. 

The final coalescence of binaries consisting of two comparable-mass black-holes, 
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with mass ratios q = M\jM% ~ (1 — 10), is expected to be one of the strongest 

astrophysical sources of energy the form of gravitational radiation. As the emis- 
sion of gravitational waves removes energy and momentum from the binary, the 
orbits shrink and the black holes eventually merge together into a single black 
hole, producing an intense burst of radiation. With the advent of ground-based 
gravitational wave detectors such as LIGO and VIRGO (which will detect merg- 
ers of black holes with masses in the range ~ (10 — 100)Mq), and with planning 
underway for the space-based LISA (which will observe mergers of massive black 
holes, with masses ~ (10 4 — 10 6 )Mq), knowledge of black-hole binary gravita- 
tional waveforms is urgent. 
1.0.3 Black-Hole Binary Coalescence Black- hole binary coalescence 
can be thought of as proceeding in three stages. During the inspiral the holes 
have wide enough separations that they can be treated as point particles. In this 
stage, the orbital period is much shorter than the timescale on which the orbital 
parameters change, and the holes spiral together on quasi-circular orbits. When 
the holes get so close together that they can no longer be approximated as point 
particles, they enter the merger phase. In this strong-field, dynamical regime, the 
black holes plunge together and merge into a single, highly distorted black hole, 
surrounded by a common horizon. This remnant black hole then "rings down," 
shedding its nonaxisymmetric modes through gravitational wave emission and 
settling down into a quiescent rotating black hole. 

1.0.4 Gravitational Waves from the Coalescence Knowledge of the 
gravitational waveforms from these three stages of black hole coalescence is im- 
portant for gravitational wave detection and data analysis, as well as astrophys- 
ical applications. The inspiral can be calculated analytically using the post- 
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Newtonian (PN) approximation, which is an expansion of the full equations of 

General Relativity in powers of e ~ v 2 /c 2 ~ GM/(Rc 2 ), where v is the charac- 
teristic velocity of the source, M is its mass, and R is its characteristic size. The 
inspiral gravitational waveform is a "chirp" , which is a sinusoid increasing in both 
frequency and amplitude. The ringdown can also be calculated analytically us- 
ing techniques of black-hole perturbation theory, and the resulting gravitational 
waveforms are damped sinusoids. However, the merger stage can only be under- 
stood using full numerical simulations of the Einstein equations in 3-D, and the 
resulting gravitational waveforms were unknown - until recently. 

The final merger of comparable-mass black holes is a powerful source of grav- 
itational radiation. The gravitational waveforms emitted by a black-hole binary 
with mass ratio q scale with the total mass M of the binary. Setting c = 1 and 
G = 1, we can express both length and time scales in terms of the total mass: 
1M ~ 5 x 1O~ 6 (M/M )sec ~ 1.5(M/M )km. The strong-field merger will pro- 
duce a burst of gravitational radiation lasting ~ lOOAf and having a luminosity 
~ 1O 23 L , which is greater than the combined luminosities of all the stars in the 
visible universe. For stellar black- hole binaries (M ~ 1OM ) this luminosity will 
last for ~ 5 ms, and for massive black-hole binaries (M ~ 1O 6 M ), for ~ 10 mm. 

1.0.5 Calculating the Merger Understanding the final merger of a 
black- hole binary requires solving the full Einstein equations using the methods of 
numerical relativity. This endeavor has proved to have many challenging aspects, 
ranging from providing astrophysically relevant initial data to understanding the 
structure of the Einstein equations, and the solution eluded researchers for many 
years. Recently, a series of dramatic breakthroughs has ignited the field, mak- 
ing robust, stable, and accurate simulations of binary mergers possible for the 
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first time. These models are opening our understanding of strong-field dynamics, 

impacting gravitational wave searches and other areas of astrophysics. In this 

article we provide an overview of these exciting developments, highlighting the 

key physical results that are emerging. 

2 Steps toward a robust black-hole binary model 

Computing the strong-field merger of two comparable-mass black holes has a 
long history, with the first attempt dating back more than 40 years. Overall, 
progress was generally slow and incremental, requiring the interplay among new 
ideas in black hole modeling, mathematical investigations into the structure of 
the Einstein equations, the development of effective gauge conditions, advances 
in computational techniques, and the availability of high performance computing 
resources. Here we provide a brief review of key developments in this arena. 

2.1 Overview of Numerical Relativity Issues 

Solving Einstein's equations on a computer typically requires slicing 4-D space- 
time into a stack of 3-D spacelike hypersurfaces, each labeled by time t dUEJED. 
The Einstein equations then divide into two sets: constraint and evolution equa- 
tions. The constraints are a set of time-independent elliptical equations that 
must hold on each slice. In particular, the constraints are solved first to obtain 
valid initial data for a black-hole binary simulation. This data is then propagated 
forward in time using the evolution equations. 

2.1.1 Gauge Freedom General Relativity has four spacetime coordinate 
degrees of freedom, which are associated with four freely-specifiable coordinate or 
gauge functions that govern the future development of both time and the spatial 
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coordinates. These are generally taken to be the lapse function a, which gives 

the lapse of proper time (adt) between neighboring slices, and the shift vector 

f3 l , which determines how the spatial coordinates develop from one slice to the 

next. Appropriate choices for the lapse and shift have proven to be crucial for 

successful black hole evolutions (JH[5|). 

2.1.2 Formulation of the Evolution Equations For many years, a 
primary challenge of numerical relativity was simply to decide which equations 
to solve ([3]). There is no unique formulation of the Einstein equations; rather 
there are many choices regarding, for example, which variables to use, whether 
to write the equations as first-order or second-order in time, and which coordinate 
conditions to impose. These choices are not arbitrary because some formulations 
turn out to be more "numerically friendly" than others. That is, although they 
are analytically equivalent when constraints and auxiliary variable definitions 
are assumed exact, some formulations may be unstable in practice. For exam- 
ple, the evolution equations may admit rapidly growing solutions which violate 
the constraints. In this case, although the initial data may contain only tiny 
errors, the subsequent evolution may produce violations of Einstein's equations 
which "blow-up". There can also be pathological numerical solutions that don't 
represent solutions to the analytic evolution equations at all, but are supported 
by the discrete numerical grid structure and depend on the details of the finite- 
differencing operators. 

2.1.3 Numerical Methods Once a formulation has been chosen, the Ein- 
stein equations are solved numerically for various field variables on a grid of 
discrete points that represents the spacetime domain of interest. There are two 
general approaches for dealing with the spatial derivatives that appear in these 
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equations. Finite differencing interpolates the derivative at a given point from the 

surrounding points according to a Taylor expansion in the grid spacing. Spectral 
methods assume a solution in the form of a summation of orthogonal functions; 
once the coefficients are obtained numerically, the derivatives can be found an- 
alytically. In both cases, time integration can be handled in a number of ways, 
most commonly by a Runge-Kutta algorithm. To date, most numerical relativity 
solutions have been carried out using finite differencing, with most current work 
employing a Cartesian grid in three spatial dimensions, although results from 
evolutions using spectral methods are now becoming available. 

2.1.4 initial conditions FOR black holes To model astrophysical bina- 
ries, initial data must be generated for two black holes moving on quasicircular 
inspiralling orbits and having a mass-ratio q and some configuration of spins. 
One usually conceptualizes the initial state of the system with a particle-like de- 
scription of the masses, positions, momenta, and spins of each black hole. Since 
general relativity is a field theory, such a description can only be seen as a first 
step. The model requires a full description for the initial field configuration which 
satisfies General Relativity's four initial value constraints, and which somehow 
corresponds to the system we have described in these particle-like terms. 

2.1.5 Evolving the Binary This system is evolved for several orbits and 
then through plunge, merger, and ringdown, for a total duration on the order of 
several hundred M or more. To obtain the gravitational waveforms, the radiation 
must be extracted from the simulation at large enough distances from the source 
to be in the "wave zone". Since the length scales of the black holes are ~ M 
and the wavelengths of the emitted radiation are ~ (10 — 100)M, it is clear that 
some sort of variable resolution such as adaptive mesh refinement must be used 
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to handle the very large computational domains needed. 
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2.2 The Lazarus Approach: Hybrid Black Hole Merger Wave- 
forms 

By the late 1990's numerical relativity had developed to the point that brief simu- 
lations of black-hole binaries in three spatial dimensions plus time were possible. 
These techniques were sufficient for evolving promptly merging "grazing colli- 
sions" of black holes (JHIED- However, the simulations were not indefinitely stable, 
but rather would typically crash after ~ (10 — 30) M, well before any significant 
portion of a binary orbit could be evolved. 

2.2.1 Hybrid simulations In this arena, the Lazarus method emerged as 
a novel approach to obtaining black-hole binary waveforms, combining short nu- 
merical relativity simulations with black hole perturbation methods (|8|). Since 
numerical relativity codes were then only able to evolve for a brief period, and 
since perturbation theory could approximate the late time dynamics of the dis- 
torted remnant black hole, the Lazarus Project aimed to apply numerical rela- 
tivity to evolve the strong-field approach to merger, thereby providing a hybrid 
model for a significant part of the problem. Starting from quasi-circular initial 
configurations near the start of the final merger, the black holes were evolved 
using numerical relativity for ~ 15M until just before the simulation became 
inaccurate. Then, via a complicated interface, data from late in the numerical 
simulations was interpreted as initial data describing the perturbed final black 
hole and the emerging radiation. Finally, black hole perturbation theory tech- 
niques were applied to evolve this data and derive the full waveforms. 
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2.2.2 First Results The Lazarus simulations gave the first indication of 

what might be expected for the terminal burst of radiation from coalescing black- 
hole binaries. Figure Q] shows that the dominant I = 2,m = 2 spin-(-2)-weighted 
spherical harmonic component has a brief burst of radiation smoothly joining the 
damped sinusoidal signal of the ringdown. The waveforms were remarkable for 
their simplicity, with predominantly circular polarization (in the I = 2, m = 2 
mode), and a steady evolution of polarization frequency and amplitude. Since 
these characteristics were robust under variations in the model, it was conjectured 
that the waveforms from a binary starting from a wide separation late in the 
inspiral would be well characterized by these simple features. Subsequently, the 
Lazarus method was applied to study mergers of more generic black-hole binaries 

2.3 Toward Evolving a Black-Hole Binary Orbit 

2.3.1 The BSSN system Concurrently with the Lazarus investigations, 
further progress was being made toward full numerical relativity simulations of 
black-hole binary mergers. One major milestone was the development of a con- 
formal formulation of the Einstein equations known as the Baumgarte-Shapiro- 
Shibata-Nakamura or "BSSN" system (|12p . that overcame some instability prob- 
lems associated with an earlier formulation in use at the time (fTj ll3p . In this 
approach, the set of evolution equations is written with first-order time deriva- 
tives and second-order spatial derivatives, and is strongly hyperbolic (|14p . Stable 
time evolution was accomplished using a coordinate condition that evolves the 
lapse function dynamically, causing the slices to avoid crashing into the black 
hole singularities (|15^ I7|). However the simulations were still limited to durations 
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< (30 — 40) M of stable evolution, by a failure of the spatial coordinate system 

known as "grid stretching," in which the coordinates tend to fall into the black 

holes, and by instabilities related to how the black holes were handled numerically. 

2.3.2 Controlling the spatial coordinates Eliminating grid-stretching 
required developing appropriate techniques for dynamically controlling the spa- 
tial coordinates, which are governed by the evolution of the shift vector. The 
first long-lasting evolutions of distorted black holes relied on a class of hyperbolic 
shift evolution schemes, known as T-freezing conditions, which were inspired by 
the BSSN formulation. With this approach, a single distorted black hole could 
be evolved indefinitely (e.g., for thousands of M) with reasonable accuracy (|16p . 
These studies allowed the full numerical determination of gravitational waveforms 
from black hole ringdowns as the distorted black hole settled down to a physically 
and numerically stable final state, and provided a foundation for future advances 
in black-hole binary simulations ([5|). 

2.3.3 Problems with moving black holes However, another aspect 
of spatial gauge conditions remained a critical issue for long-lasting black-hole 
binary simulations. As general relativity allows arbitrary coordinate systems, 
many groups adopted coordinate conditions which did not allow the black holes 
to move through the computational domain. This simplified the problem of 
dealing with the black hole singularities, which were handled either by excising 
the black hole interiors (within the horizons) from the computational domain (|17p . 
or by representing the black holes as punctures (|18|) . Though progress had been 
made toward implementing moving excision regions (|19p . it was computationally 
much simpler to demand that the excision region remain fixed on the grid (|20p . 
Similarly, the puncture treatment was formulated so that the black hole was fixed 
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on the grid and the singularity was factored out and handled analytically in a 

time-independent manner. However, for binary configurations in which the black 
holes physically move, the cost of keeping the black holes fixed in coordinate 
space was paid by the twisting and stretching of the dynamical fields, eventually 
leading to large computational errors. 

For inspiraling binary configurations one potential solution was to attempt to 
untwist the geometric field data by imposing comoving coordinates using a shift 
vector that is dynamically adjusted during the evolution of the binary to mini- 
mize motion of the black hole horizons (|2ip . With this approach, Briigmann et 
al. achieved a significant milestone: the first simulation of a full orbital cycle 
of a black- hole binary system prior to merger (|22p . In this work, the authors 
applied a dynamically adjusted corotating frame tracking the measured position 
of the black hole apparent horizons. This allowed a simulation which remained 
accurate for a little longer than the > 100M duration of a complete orbit. How- 
ever, problems occurred with the large characteristic speeds that resulted near 
the outer boundaries, limiting the domains of these simulations and preventing 
gravitational waveforms from being measured effectively. The coordinate control 
also required fine tuning which eventually failed, causing the simulation to crash 
as the black holes finally approached merger. 

2.4 Robust black-hole binary simulations 

2.4.1 First Orbit and Merger simulation In early 2005 Pretorius 
shocked the numerical relativity community by announcing the first complete, 
robust simulations of an equal- mass black hole merger (|23j) . After completing 
~ 1 orbit, the black holes plunged and merged to form a single distorted black 
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hole that then rings down. Pretorius extracted the gravitational waves to obtain 

the first inspiraling merger waveform directly from numerical relativity, shown in 

Fig. m 

Pretorius employed several techniques that were very different from most other 
numerical approaches to the black-hole binary problem (|24p . Rather than using 
the BSSN formulation, he applied a generalized harmonic formalism (|25p . directly 
integrating the spacetime metric with evolution equations of second-order in both 
space and time. These equations were implemented numerically using adaptive 
mesh refinement to allow high resolution around the black holes while maintaining 
a large computational domain. Pretorius utilized spatial coordinates compactified 
to draw spatial infinity into the computational domain, with a choice of gauge- 
evolution strongly tied to his evolution formalism. Following (|26p . he added 
terms to the evolution equations to specifically damp away any violations of 
the constraints. In his simulations, the black holes were excised and moved freely 
across the computational grid. Pretorius' success with such novel methods quickly 
raised questions as to whether the struggling, more widely-pursued BSSN-based 
puncture approach might be off course. 

2.4.2 Moving Punctures Later that same year, however, a new robust 
method based on the BSSN formulation was announced. The "moving punc- 
ture" method was discovered simultaneously and independently by the groups at 
the University of Texas at Brownsville (UTB) (|27|) and NASA's Goddard Space 
Flight Center (GSFC) (|28p . In this approach the black holes are represented 
as punctures, but are not constrained to remain fixed on the coordinate grid. 
Rather, they are allowed to move freely through the grid using novel coordinate 
conditions. Figure [3] shows the trajectories of the puncture black holes as com- 
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puted by the UTB group; the formation of a common horizon marks the point of 

merger. The moving puncture method eliminates the analytic representation of 

the puncture singularities in favor of an approximate numerical treatment within 

the black hole horizons. 

The UTB and GSFC groups had discovered and applied similar methods to 
the same problem: evolving an equal-mass nonspinning black-hole binary from 
near the final orbit, through merger and ringdown, and studying the gravitational 
waves. The first generation of merger waveforms from Pretorius, GSFC, and UTB 
showed the same simple burst of radiation ending in a damped-sinusoidal ring- 
down, and were qualitatively consistent with the Lazarus project results discussed 
above and shown in Fig. [TJ 

The discovery of the moving puncture method ignited the field of black-hole bi- 
nary evolutions. Since it was based on commonly used methods, most researchers 
in the field were quickly able to achieve accurate and stable evolutions using their 
existing codes, with the adoption of simple coordinate conditions (|29p . Suddenly, 
the game was on and nearly all the groups were participating. 

3 The Physics of Black-Hole Binaries 

With the advent of successful numerical evolutions of binaries that inspiral and 
merge, the numerical relativity community's focus changed to investigating the 
physics of binary mergers. This advancing frontier is enabled by continued techni- 
cal improvements in areas ranging from initial data prescriptions to more accurate 
numerical methods. 
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3.1 Merger Dynamics and Waveforms 

The astrophysical black hole mergers that are the targets of current and future 
gravitational wave detectors are expected to reach the merger stage after having 
radiated away any initial eccentricity (p?Uj) and proceeding through a long quasi- 
circular inspiral. All equal-mass, nonspinning binary merger simulations starting 
from such orbits in the inspiral should produce the same gravitational waveform, 
subject only to rescaling with the total mass of the system. For many years, 
concerns had been raised about the accuracy and realism of black-hole binary 
initial data sets, including the effects of spurious gravitational radiation and 
eccentricity (|31U32l[33|) . With numerical relativity now able to simulate the final 
merger, the next step was to run models with enough orbits before the plunge 
and merger to get complete and reproducible waveforms starting from the late 
inspiral. 

3.1.1 Equal-Mass, Nonspinning Black Holes The GSFC group pro- 
duced the first representation of the definitive waveform for the final stages of a 
merger of equal-mass, nonspinning black holes (|34|) . They carried out a series of 
four simulations with the holes starting from quasicircular initial conditions at 
increasingly larger separations. In these runs, the holes completed ~ 1.8,2.5,3.6, 
and 4.2 orbits before the formation of a common horizon. To compare the results 
of these models, they chose the moment of peak gravitational radiation amplitude 
as the fiducial time t = 0. 

The orbital dynamics of the binaries can be examined by tracking the black 
hole centers, given by the location of the punctures. The black hole trajectories 
for each run were oriented so that they superpose at the fiducial t = 0. In the 
early stages of each run, the tracks clearly showed the effects of eccentricity in 
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the initial conditions, with the amount of initial eccentricity decreasing for more 

widely separated holes. As the holes spiraled together deeper into the strong-field 

regime, the tracks locked on to a universal trajectory independent of their initial 

conditions that continued for the last orbit, plunge, and merger. 

Figure H] reveals the corresponding universal gravitational waveform. Here, the 

dominant i = 2, m = 2 quadrupolar component for each run is shown, shifted 

in time so that the peak radiation amplitude occurs at t = 0. Starting from 

t = — 50M the waveforms show nearly perfect agreement, differing from each 

other at the level of 1%. The signals for the preceding few orbits agree at the 

level of ~ 10%, except for a brief burst of spurious radiation at the start of each 

run. 

Note that the merger waveform shows a remarkably simple shape, making a 
smooth transition from the inspiral chirp to the damped sinusoid of the ring- 
down. As the merger begins, both the wave amplitude and frequency increase, 
albeit faster than in the inspiral. The amplitude reaches a peak and then de- 
creases, dropping exponentially through the ringdown. The frequency increases 
monotonically to a maximum value that remains constant during the ringdown. 

Of course, the black- hole binary merger dynamics and waveforms can be altered 
by the presence of large amounts of eccentricity ([35]) and spurious gravitational 
radiation (36) near the time of the plunge. However, the robustness of the merger 
to modest deviations from astrophysical initial conditions opened the door to 
studying many cases of interest using relatively short simulations, starting just a 
few orbits before the start of the plunge. 

3.1.2 Unequal-Mass, Nonspinning Black Holes Astrophysical black- 
hole binaries are unlikely to have exactly equal masses. Currently, numerical 
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relativists are able to evolve systems with mass ratios up to q = 10 (37). Starting 

from quasicircular orbits the simulations show that the merger phase for non- 
spinning, unequal mass black hole binaries is robust to modest deviations from 
these initial conditions and produces a generally simple waveform shape. 

An important tool for analyzing these mergers is a decomposition into spin- 
weighted spherical harmonic modes. Berti et al. (|38|) analyzed a set of unequal- 
mass nonspinning mergers with mass ratios ranging from q = 1 to q = 4. Studying 
the multipolar distribution of the radiation, they found that the sub-dominant 
modes (£ > 2) become more important, carrying a larger fraction of the energy, 
as q increases. Specifically, for q > 2, the t = 3 mode typically carries ~ 10% of 
the radiated energy. Also, as expected from symmetry considerations, the odd-m 
modes are suppressed in the equal-mass limit. 

Baker et al. (|39p carried out a complementary study of the radiation from 
nonspinning mergers with mass ratios in the range 1 < q < 6. The multipolar 
decomposition clearly shows that the hallmark simplicity of the waveform persists 
for q > 1 and extends to each of the spherical harmonic components I > 2; this 
property has recently been shown to extend to the q = 10 case (|37p . In the 
full mode-summed waveform, this simple shape is also seen when viewing along 
the system's orbital axis, where the quadrupole mode dominates. A somewhat 
more complex appearance arises by viewing the system off- axis, where higher 
multipoles contribute more strongly to the waveform at various angles. 

Throughout the entire coalescence, each of the spherical harmonic waveform 
components is circularly polarized, with steadily varying phase and amplitude 
(|39p . For each mass ratio q, the rotational phase (and frequency) of the different 
(£, m) components are the same. During the inspiral this is expected, since the 
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waveform phase is equal to the rotational phase multiplied by the mode number 

m. However, for the £ = m modes this relationship also holds throughout the 

merger and into the ringdown. 

These properties suggest a simple conceptual interpretation in which the radi- 
ation is generated by an "implicit rotating source." In this picture, each (£, m) 
mode is generated separately by the (£, m) moment of some implicit source. The 
fixed relationship for the £ = m modes implies that these components of the 
source rotate rigidly through the entire coalescence. The £ 7^ m components are 
less rigid and peel away from the main £ = m trend during the merger (|39p . 
3.1.3 Mergers of Black Holes with Spin The mass ratio q is a one- 
dimensional cut into the parameter space of black-hole binaries. The remaining 
parameter space is dominated by the spin angular momentum of each hole. As 
spin is a vector quantity, this adds six more dimensions. Simulations of spinning 
black holes first focused on systems whose spins were expected, on the basis of PN 
arguments, to have the least effect on the orbital motion. Binaries with aligned 
and anti-aligned spins are relatively easy to treat as there is no precession of 
the orbital plane or the individual spins. When one or both holes has a spin 
not parallel to the orbital axis, spin-orbit and spin-spin interactions will cause 
precession that can complicate the evolution and the resulting waveforms (|40p . 

The first merger evolutions of spinning black holes were carried out by Cam- 
panelli et al. (|4ip . They simulated the mergers of two equal-mass highly spinning 
black holes with {a/M)\^ = .757, and both spin vectors either aligned or anti- 
aligned with the orbital angular momentum. Here, a < M is the magnitude 
of the black hole spin angular momentum per unit mass. They also evolved a 
nonspinning equal-mass binary for comparison. All three binaries had the same 
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initial orbital angular frequency corresponding to an orbital period ~ 125M, and 

merged to form a rotating remnant black hole with {a/M)^ n& \ < 1. However, 
the aligned system took noticeably more orbits to merge than the others. This 
behavior is caused by the spin-orbit interaction, which produces an effective force 
between the black holes, either an attractive or repulsive for the anti- aligned or 
aligned cases, respectively. All three binaries generate remarkably similar grav- 
itational waveforms having a simple shape, with the aligned case displaying a 
longer wavetrain and the anti- aligned case a shorter one. 

The first fully numerical evolutions of strongly precessing spinning systems 
were carried out by Campanelli et al. (|421 l43 p . who observed both significant 
precession of the orbits and a "spin flip," in which the final post-merger black 
hole spins in the direction opposite to the two initial black holes. More recent 
work is beginning to probe the effects of spin precession on waveforms (|44p . 

Studies of this most general of (non-eccentric) parameter sets are still in the 
early stages. No systematic study of waveform shapes and polarizations has been 
carried out yet. A preliminary study of the multipolar structure of gravitational 
waves from several classes of binaries with equal spins has been performed by 
Berti et al. (|45p . They considered an equal-mass case with aligned spins; sev- 
eral q = 4 binaries with antialigned (down-down) spins; and three unequal- mass 
binaries with spins initially in the orbital plane and pointing in opposite direc- 
tions. Examining the distribution of gravitational- wave energy emitted by various 
modes, they find that, as in the case of nonspinning mergers, odd-l multipoles are 
suppressed for q = 1 and that, as q increases, more energy is radiated in higher-^ 
multipoles. 
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3.2 The Spin of the Final Black Hole 

Because the state of the final black hole formed from the coalescence depends 
primarily on termination of the inspiral and the late burst of radiation in the 
merger, it can be accurately probed with relatively short simulations of just a 
few orbits. 

The merger of two equal-mass nonspinning black holes produces a final black 
hole with a moderately high spin, (a/M) final ~ 0.69 ^^B^^. Since the 
black holes each start out with no spin, and any tidal spin- up is negligible (|47p . 
the spin of this final black hole arises from the orbital angular momentum of the 
original binary. Simulations show that, for this simplest black hole merger, the 
final spin is "universal," or independent of the initial black hole separation, for 
modest deviations from quasicircular initial conditions. 

For mergers of nonspinning black holes with unequal masses, the spin of the 
final black hole decreases as the mass ratio q increases. Simulations with mass 
ratios up to q = 10 show that the final spin parameter scales as (a/A/)fi na i ~ 
q/{l + q) 2 , where (a/M) flna i 0.48 for q = 4 and (a/M) fina i « 0.26 for q = 10 

fl3SGEU3B- 

The effects of spin-orbit and spin-spin coupling can become important in de- 
termining the spin of the final black hole. In the simplest cases, the spin vectors 
are parallel to the orbital angular momentum. Depending on the mass ratio and 
the black hole spin, the merger can result in a final black hole with a larger (spun 
up) or a smaller (spun down) spin than either of the progenitors (|10y 4 Hl4"8]) . In 
certain cases, the merger can lead to a spin flip with the final black hole spinning 
in a direction opposite to the spins of the initial holes; in particular, it is possible 
to produce a final nonspinning black hole, (a/M)fl na i = 0, from the merger of two 
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spinning holes (|49|J45p , More general black-hole binaries, with misaligned spins, 

are further complicated with the effects of precession (42.50,51). 

Several attempts have been made to produce expressions for the final spin 

vector using analytic techniques or by fitting to results from numerical simulations 

(|49U 52 1[53|) ; see ()54|) for a review. Of particular interest is the question of whether 

a black hole merger can produce a maximally spinning hole ((o/M)fl Iia i = 1) or 

indeed exceed the Kerr limit ((a/M)fi na i > 1). Current research suggests that 

this is not possible, even for mergers occurring from hyperbolic encounters (|55p 

and mergers of highly boosted black holes ()56|) . 

3.3 Recoil Kicks from Gravitational Radiation 

A notable phenomenon arising from asymmetric binary systems is the merger 
recoil or kick - a net movement of the end-state black hole from the system's 
center of mass, caused by the anisotropic emission of gravitational radiation dur- 
ing the coalescence. Fitchett f)57[) produced a useful formula for the kick velocity, 
which has significant applications in astrophysics, using a quasi- Newtonian ap- 
proximation. Several authors also calculated the recoil analytically using PN 
approximations (58). However, since the dominant part of the effect builds up in 
the strong-field regime close to merger, full numerical relativity simulations are 
needed for accurate calculations of the kick. 
3.3.1 Kicks from nonspinning black hole mergers In 2006, recoil 
from a fully numerical binary merger was demonstrated for the first time by 
Herrmann et al. (|59p . for plunging black- hole binaries with mass ratios as large as 
q ~ 3.1. This was followed soon afterwards by a full orbit and plunge simulation 
by Baker et al. (|6Q|) . who found a recoil of between 86 and 97 kms~ 1 for a mass 
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ratio q = 1.5. 

A more systematic study of recoil from mergers of unequal-mass binaries was 
produced by Gonzalez et al. in 2007 ([6T|) . who studied the merger of binary 
systems with mass ratios between q = 1 and q = 4. Figure [5] shows the resulting 
range of recoil speeds, together with several earlier numerical and analytical es- 
timates. The authors synthesized these into a single recoil formula, a nonlinear 
correction to the Fitchett formula, yielding a maximum recoil of 175 kms -1 for a 
mass ratio of q ~ 3. This has recently been tested for the more extreme q = 10 
mass ratio, with general agreement (|37p . 

3.3.2 Spinning black hole mergers and superkicks Calculations of 
recoil in the much larger parameter space of spinning binaries began with the 
non-precessing cases of holes with spins aligned (or anti-aligned) with the orbital 
angular momentum. Several studies of this region of parameter space (|62ll48p have 
revealed that the kick velocity has a quadratic dependence on initial spins, with 
a maximum kick of 448 kms~ 1 for extremal Kerr holes, (a/M)i2 = 1. For these 
anti/aligned well as for nonspinning unequal-mass black-hole mergers, 

the direction of the kick velocity is always in the orbital plane. 

Meanwhile attention turned to more general black hole spins. Campanelli et al. 
(|43p speculated from PN arguments that optimal spin configurations could give 
rise to huge "superkicks," with velocities > 1000 kms -1 , out of the initial orbital 
plane. The first such superkick — around 2500 kms -1 — was soon observed by 
Gonzalez et al. (|63|) . Tichy et al. (|50|) have argued that superkicks arise in 
mergers with general spin orientations, while greater insight into the mechanism 
of these kicks has been developed (|64U65p . 

With such a large parameter space to cover, it seems useful to try to construct a 
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general formula that will describe kicks from arbitrarily spinning binaries. Baker 

et al. (|66p used new and existing aligned-spin results to produce a single unifying 

formula for in-plane kicks; Campanelli et al. P31l67|) proposed an extension to 

this model, with scaling of the superkick out-of-plane contribution motivated by 

PN theory. The leading dependence in this formula on the angles between spins 

and linear momenta of the pre-merger holes has strong support from numerical 

simulations (67,64,68); however, the dominant scaling with mass ratio is still in 

dispute (|69l[70|) . 

3.4 Longer Waveforms: Modeling the Late Inspiral 

Many key features of black-hole binary interactions can be modeled usefully with 
only a small handful of binary orbits before merger. Quantities such as radiated 
energy and momentum are bulked near the merger and have been shown to be 
robust to the addition of one or two extra orbits (see Section [3.1.ip . However, 
optimal observational studies of gravitational waveforms ultimately require theo- 
retical predictions for the full waveform, starting at large separations during the 
inspiral. 

Before the breakthroughs in numerical relativity, most information about the 
dynamics of compact binaries came from PN theory. This approach supplied the 
particle trajectories, energy flux, and - most importantly for detector scientists 
- the gravitational waveforms themselves. However PN theory fails before the 
system merges, so the waveforms are necessarily incomplete. With the advent of 
numerical simulations encompassing many orbits, scientists finally have a way to 
develop complete information about the waveforms. This requires longer simula- 
tions which reveal the last part of the binary inspiral, and allow overlap with PN 
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waveform predictions. 

3.4.1 Low-eccentricity initial data Simulations for these studies re- 
quire careful attention to the initial configuration of the black holes. For a sig- 
nificant population of astrophysical binaries, it can be expected that the orbits 
have circularized through gravitational-wave emission over many orbits prior to 
merger. Numerical simulations try to mirror this expectation by selecting initial 
momenta consistent with near-zero eccentricity for the relatively small separa- 
tions at which a full numerical relativity simulation becomes feasible. 

For equal-mass, nonspinning binaries, two basic approaches have proved effec- 
tive in reducing spurious eccentricity. The more direct method is to model the 
observed eccentricity, and then adjust the momenta using a Newton-like step to 
zero it out; this approach has led to extremely low eccentricities (|71t l46 p . An 
alternative approach, adopted by Husa et al. (|72p . is to use the PN Hamilto- 
nian equations of motion to model the early evolution of the particle trajectories, 
starting from large (~ 50 — lOOAf ) separations; the emission of radiation during 
this process naturally circularizes the orbit, and low-eccentricity momenta can be 
read off at the desired separation. This approach has recently been tested with 
spinning binaries as well (|44p . 

The first long waveform was produced by the GSFC group (|73|, l74 p for equal- 
mass, nonspinning black holes starting ~ 7 orbits or ~ 14 gravitational wave 
cycles before merger. Baker et al. estimated their numerical errors in waveform 
phase from (|74p as a function of frequency, finding that above a certain frequency 
numerical errors are smaller than internal errors in the PN sequence. In addition, 
the numerical and PN waveform phases agree to within one radian of phase 
drift for a little over ten gravitational wave cycles preceding the last orbit before 
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merger, comparable to numerical error estimates. Hannam et al. (75) improved 

on this work, with phase and amplitude comparison between their low-eccentricity 

higher-resolution evolutions and PN waveforms. 

3.4.2 Evolutions with Spectral Techniques Numerical simulation 
codes based on pseudospectral differencing techniques are particularly well-suited 
to long waveform studies. The Caltech-Cornell group adopted the constraint- 
damped generalized harmonic formalism used by Pretorius, an important com- 
ponent in developing a stable spectral evolution code ()76|) . Also like Pretorius, 
their code handles black holes by excising the black hole interiors. Their spectral 
code also employs a numerical grid that tracks the black holes explicitly. Though 
there can be difficulties with the changing geometry as black holes merge, this 
approach allows efficient study of the long-lasting inspiral waveforms. 

Recently this approach has provided, several of the longest and most accurate 
black-hole binary evolutions . (|46 II 77j) In particular, the Caltech-Cornell group have 
used their spectral code to simulate an equal-mass nonspinning binary starting 
16 orbits and 32 gravitational wave cycles before merger: see Fig. [6] (|46p These 
much longer waveforms have been used to validate several competing PN models 

3.4.3 Comparing results With so many pre-merger waveform cycles now 
available for the equal-mass case, the results can be cross-checked by compar- 
ing the "complete" waveform - inspiral, merger and ringdown - between groups. 
This is important to verify expectations that differences in methodology and 
residual numerical effects, such as unwanted eccentricity, are unimportant. A re- 
cent effort, dubbed the "Samurai project," analyzed long waveforms from several 
groups in the light of detectability criteria for the LIGO and Virgo ground-based 
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gravitational-wave detectors. They found that the available numerical relativ- 
ity waveforms are indistinguishable in these detectors for signal-to-noise ratios 
(SNRs) < 25 CZ5J- 

3.4.4 Long Waveforms For More General Black Holes Hannam 
et al. have investigated the properties of highly spinning orbit-aligned black 
holes, comparing their phase and amplitude with PN predictions over the last 
ten waveform cycles (|79p. In the best cases, they find phase agreement within 2.5 
rad for 3.5PN, and amplitude agreement to within around 12% with restricted 
PN. While not at the same level as their nonspinning results (|75j). this may be 
attributable to the lower orders of accuracy available for spinning black holes in 
PN theory at the present. The Caltech-Cornell group has recently conducted 
long-lasting simulations of black holes with spins aligned and anti-aligned with 
the orbital angular momentum, calibrating a tunable PN waveform model to 
match the results (|77 H80j) ; see Sec. 15.11 

Generic mergers involving non-aligned, and thus precessing, spins adds four 
new degrees of freedom to the problem. A systematic understanding of wave- 
forms generated by generic mergers will require considerably more study. Recent 
simulations have begun to explore generic examples (|44p . 

4 Applications in Astrophysics 

We have described a sampling of the explosion of numerical relativity studies 
revealing some of the details of black-hole binary physics as implied by General 
Relativity. While more remains to be learned, this new understanding is already 
making important contributions in planning and interpreting astrophysical black 
hole observations where Einstein's theory is applied as the standard model of 
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gravitational physics. 
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4.1 Waveforms for Gravitational- Wave Observations 

Experimental gravitational-wave detectors were first developed more than 40 
years ago. Although advances in design have increased their sensitivity by many 
orders of magnitude, the extreme weakness of expected astrophysical signals 
(strain amplitudes 5L/L ~ 10~ 21 ) means the observational challenge is still huge. 

The output of any gravitational wave detector will be a data stream that must 
be combed through to find real signals. This search requires accurate "template" 
waveforms representing our best picture of the radiation from expected sources; 
these templates can then be compared with the data stream through matched 
filtering. We refer the reader to a review on gravitational-wave astronomy for an 
overview of these techniques (j8Tj) . 

Crucially, the most important sources of gravitational radiation are expected 
to include the mergers of black-hole binaries. Before the advent of numerical rela- 
tivity simulations of black hole mergers, the only test waveforms available for use 
in data analysis studies were based on PN theory. These waveforms were essen- 
tially only inspiral chirps and did not include the strong-field merger. The new, 
richer, information now available from numerical relativity has revolutionized the 
data analysis picture in several ways. 

4.1.1 Detecting black-hole mergers The availability of the plunge- 
merger-ringdown portion of the signal can greatly increase the SNR in the detec- 
tor. Armed with a long numerical merger waveform of acceptable accuracy, we 
can extend it backwards to cover an arbitrarily long inspiral by matching to a 
PN waveform. Such a "hybrid" waveform was first produced by the GSFC group 
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for the equal-mass, nonspinning case (|74p . Using this hybrid, we can investigate 

the total achievable SNR, and the related distance reach, for current and future 
detectors. Figure [7] demonstrates the gain in SNR from including the merger por- 
tion of the waveform for the ground-based Advanced LIGO detector. Contours 
of SNR as function of redshift z and total binary mass M for the LISA detector 
are shown in Figure [HJ 

Full numerical waveforms, and the longer hybrid waveforms generated from 
them, can also be used to improve existing data-analysis techniques and template 
sets. Since previously developed gravitational wave data analysis algorithms were 
not based on knowledge of the merger waveforms, an obvious first step is to test 
how successfully these techniques detect the merger signals predicted by numeri- 
cal simulations. In 2009, the NINJA project (|8"2"j) used direct injection of a broad 
range of short and long numerical merger waveforms into mock LIGO and Virgo 
data streams for this purpose. The result was the most realistic testing ground 
to date for disparate data- analysis methods, including full and partial waveform 
template matching, as well as unmodeled burst searches. Further studies of de- 
tection algorithms with numerical relativity waveforms are continuing. 
4.1.2 Measuring black-hole binary parameters The SNR is only a 
crude guide to the specific detector response to gravitational waveforms, however. 
The merger portion of the waveform, though short in duration, may contain 
important new information not present in the inspiral. As a simple example, 
we expect the time of merger itself to be well localized with the full waveform, 
whereas it is not well-defined in the inspiral-only signal. 

More generally, the different (I, m) modes of a binary's full inspiral-merger- 
ringdown waveform scale differently with the time to merger. Modes that were 



Final Merger of Black-Hole Binaries 29 
not significant during inspiral suddenly become more prominent in the merger, 

and the detailed information they carry becomes available to the observer (|38p39p 
After detection of gravitational waves from distant sources, we are most inter- 
ested in identifying the physical parameters of the sources. Each of the seven 
intrinsic parameters of a black-hole binary (mass ratio and spin vectors) will, in 
general, be imprinted on the gravitational-wave signal, along with some extrin- 
sic parameters such as sky position and distance to the source. For sufficiently 
strong SNRs, expected for massive binary mergers seen by LISA, we can expect 
to be able to extract some of these parameters at high precision. While these pa- 
rameters can be partially disentangled using inspiral-only template information, 
it has recently been found that the full merger waveform can help break parame- 
ter degeneracies and hence drive down uncertainties in several important physical 
parameters. Good localization of the source on the sky is especially important for 
the development of multi-messenger astronomy. The recent parameter estimation 
studies of non-spinning mergers that include the merger waveforms indicate that 
LISA will be able to locate sky positions within a few arcminutes for binaries 
with ~ 10 6 Mq at cosmological distances (redshift z = l) (|83p . 

4.2 Consequences of Merger Recoil 

As previously discussed, asymmetries in a black-hole binary system due to un- 
equal masses and/or spins result in the anisotropic emission of gravitational ra- 
diation, ultimately imparting a recoil to the merged remnant black hole. Numer- 
ically it has been found that, for certain configurations of black hole spins, the 
recoil velocity can exceed the escape velocity of many galaxies. Since it is impor- 
tant to determine the theoretical probability that a massive black hole might be 
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ejected from its host galaxy, there have been some preliminary calculations based 

on simple distributions of spins and mass ratios, and analytic fits of numerical 

results giving the dependence of the recoil on mass and spin (|84t l69 p . Although 

a consensus on the exact probability of galactic ejection has yet to be achieved, 

there is general agreement that it is non-negligible. 

Such rogue black holes may have already been observed in the form of two 
rapidly moving quasars (|85p . It has been speculated that these particular quasars 
originated from the mergers of massive black holes during the coalescence of their 
host galaxies. However a previous study of quasar data found no indications of 
such recoil events, suggesting that they are rare (86). 

The effect of recoil may also be observed less directly in its effect on the growth 
rates of black holes. Those ejected into the sparse intergalactic medium are less 
likely to encounter and merge with other black holes (|87|) , Even the growth of 
recoiled black holes remaining within their host galaxies may be affected, as the 
motion of the black hole can modify the rate at which it accretes matter (j88[) . 

4.3 Mass and Spin Evolution 

The expected distribution of masses and spins of astrophysical black holes is an- 
other topic of considerable interest in astrophysics. In our current understanding, 
black holes grow from smaller "seeds" early in the history of the universe through 
a combination of mergers and the accretion of gas (|89p . In general, most of the 
mass growth is believed to come from accretion, with mergers providing a modest 
increase. 

However, the gravitational radiation emitted during black hole coalescence car- 
ries away energy, reducing the overall system mass by roughly several percent. 
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The bulk of this loss happens quickly, in the final plunge and merger stage of 

the coalescence. With this rapid mass loss, nearby matter in an accretion disk 

around the remnant black hole might react to the accompanying sudden change in 

gravitational potential and produce a visible change in its electromagnetic profile 

- a possible electromagnetic counterpart to the burst of gravitational radiation 

dsnnHU). 

For the final spin the situation is considerably more interesting, and may an- 
swer questions about double-jet "X-shaped" radio sources (92). We have seen in 
Sec. 13.21 that the final spin from black holes merging in vacuum depends on the 
magnitudes of the initial spins and their orientations relative to the orbital axis. 
The RIT group (j53]) has recently used their extended mass and spin formulas 
in a spin-evolution study, obtaining asymptotic spin distributions for BH merger 
remnants assuming no accretion. However, interaction of the merging holes with 
surrounding gas may serve to align the binary spins before merger, changing the 
picture somewhat (|93p . In general, the effects of binary merger and accretion 
have to be studied together for a coherent picture to emerge (|94p . 

5 Outlook 

5.1 Complete Analytic Waveform Models 

Gravitational wave observatories may be sensitive to hundreds or thousands of 
wave cycles. Analysis of the observed data requires comparisons with model sig- 
nals representing the full variety of possible mergers. Computing so many orbits 
of pre-merger evolution using numerical relativity would be computationally very 
expensive. 

The PN approximation provides accurate representations of the dynamics and 



32 Final Merger of Black-Hole Binaries 

waveforms for the long-lasting inspiral portion of the coalescence during which 

the black holes remain fairly well-separated and their velocities remain relatively 
low. The most valuable waveform models for gravitational wave data analysis 
must combine the efficiency of the PN approach, while accurately representing 
the final merger portion of the which is only understood by numerical simulations. 
This requires a means of analytically encoding the merger signals. 

Several approaches are currently being explored for constructing these complete 
signal models. Some are based on the analytic "effective-one-body" model of 
binary coalescence (|95 U 96 |l9T|) , Figure [9] shows a comparison between an effective- 
one-body waveform for the merger of a nonspinning black hole binary with mass 
ratio q = 4 with a numerically simulated waveform ()96|) . For specific cases, such 
models have been tuned to high accuracy to agree to agree with numerical results 
(|98p . Another approach, which models the phenomenological shape of hybrid 
waveforms in frequency space (99), has been extended for the dominant waveform 
modes from spinning, but non-precessing, mergers (jlOOp . Further development 
of these models, and the production of a family of simulated waveform spanning 
the parameter space, is a current focus of broad-based research collaborations. 

5.2 Improved Numerical Methods 

5.2.1 Initial Data The initial data models currently used to begin numer- 
ical simulations do not perfectly represent the intended astrophysical systems. 
Research continues to improve these models. For example, although some as- 
trophysical black holes are expected to have near-extremal spins, common initial 
data cannot represent holes with (a/M)j > 0.93 (|101|.I102]) ; novel methods are 
being developed to go beyond this limit (jlOip . Also, current initial data mod- 
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els generally do not contain the physically appropriate radiation content for an 

inspiraling binary (|32 |,I103[) . but rather harbor spurious radiation that is not as- 
trophvsical (|104p . Efforts are underway to include initial radiation content more 
consistent with PN predictions (j!05|) . In addition, in the case of moving punc- 
tures, the initial coordinates typically undergo a rapid transition at the start of 
the evolution as the black holes relax into a more stable solution. While the 
resulting transient pulse of "gauge radiation" does not alter the physics, it does 
contain fairly high frequencies that can be challenging to resolve and it has moti- 
vated construction of analytic initial coordinates that are closer to the numerically 
evolved coordinates (|106p . 

5.2.2 Evolution More efficient and accurate methods of numerical evo- 
lution would make simulations of many binary orbits, particularly with large 
mass ratios, more computationally practical. For representing spatial derivatives 
on a computational grid, the highest accuracy of finite differencing stencil yet 
achieved in the context of numerical relativity is 8th order in the grid spacing 
(although the accuracy of interpolation between refinement boundaries is not 
yet commensurate) (|107p . Spectral methods are generally more accurate but 
less robust than finite differencing, requiring fine tuning for generic black hole 
binaries i! ION;;: they are also unlikely to handle shocks in accreting matter (see 
Section f5.3|) . Alternatives to both finite differencing and spectral methods, such 
as finite element methods, are also being investigated (|109p . Meanwhile more 
efficient time-integration techniques allowing larger step sizes are being explored 
(HID]). 

5.2.3 Wave Extraction Because the physical domain of a typical simu- 
lation is finite, gravitational radiation is usually extracted on a sphere of finite 
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radius rather than at infinity. If multiple extraction surfaces are employed in a re- 
gion of sufficiently high spatial resolution, then the radiation can be extrapolated 
to spatial infinity. A more accurate method known as "Cauchy-characteristic 
extrapolation" extracts the gravitational wave data at a finite radius and then 
evolves it along null geodesies to future null infinity; this method is currently 
under development (|llip . 

5.3 Including Gas and Magnetic Fields 

In addition to being gravitationally "loud" , black hole mergers may also be elec- 
tromagnetically visible. Massive black holes at the centers of galaxies are typ- 
ically surrounded by gaseous accretion disks and magnetic fields. When the 
black holes merge, the dynamics of the gas and magnetic fields may produce 
electromagnetic signals, counterparts to the emitted gravitational radiation. For 
example, the inspiral may "twist" the fields, resulting in characteristic electro- 
magnetic radiation (|1121lll3"j) as well as heating of a surrounding accretion disk 
(|114p . In addition, the violent merger dynamics may induce shock waves in ac- 
creting matter, in turn generating electromagnetic radiation. The recoil of the 
merged remnant, in particular, may have such an effect on the accretion disk 

jnainsiinziisBinsiinaisiiiiaij . 

5.3.1 Multimessenger Astronomy Detection of electromagnetic coun- 
terparts of gravitational waves would be of great scientific value. Current models 
of the complex merger physics (e.g. (|12ip ) could be directly tested. Einstein- 
Maxwell theory, the coupling of gravitational and electromagnetic fields on macro- 
scopic scales, could be verified. In particular, equality of the speed of grav- 
ity with the speed of light could be confirmed (|118pil2|) . Astronomy would 
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also benefit enormously, as the location and characteristics of gravitational wave 

sources could be corroborated. In addition, electromagnetically visible mergers 
could serve as "standard candles", beacons by which to measure the accelerat- 
ing expansion of the universe, while simultaneously playing the role of "standard 
sirens" emitting gravitational radiation. Cross-correlating these signals could re- 
sult in measurement of the cosmological "dark energy" to unprecedented accuracy 

C22iinsirainsiiT2u . 

5.3.2 MODELING MATTER There has been some preliminary work on the 
dynamic effects of the spacetime of a coalescing binary on surrounding matter. 
Modeling the accretion disk as geodesic particles, large collision energies were 
found as the binary merged (|125|) . And hydrodynamically simulating a gas cloud 
around around the binary, luminosity due to shocks has been calculated (|126[) . 
The generation of electromagnetic radiation by more direct means has also been 
simulated, via the twisting of a magnetic field anchored in the accretion disk, 
as it is frame-dragged by the binary (|127U113llll2p . Future efforts will employ 
magnetohydrodynamic methods, where challenges include adequately resolving 
shocks and turbulence in a dynamic spacetime, and accurately representing the 
divergence-free magnetic field on an adaptively refined grid. 

6 Conclusion 

We hope to have conveyed some of the excitement of recent progress in under- 
standing black- hole binary physics, and the applications of this knowledge in 
astrophysics. These advances are the result of sustained efforts by a broad scien- 
tific community over many years. In a brief review, it is impossible to adequately 
represent all of the excellent work that has contributed to the current state of 
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knowledge. We have only been able to provide a few of the highlights as seen 

through the lens of our particular perspective. 

We encourage interested readers to pursue the subject further. Other resources 

are available on topics including: numerical relativity techniques (|128|, [3"]). the 

breakthroughs in black hole merger simulations (j!29j) . black hole simulations for 

gravitational wave data analysis (|130j ), and gravitational wave science generally 
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t/M 



Figure 1: Waveforms from Lazarus models of equal-mass, nonspinning black hole 
mergers. This figure shows the I = 2,m = 2 component of Re(r^4), corre- 
sponding to one of the two polarization states of the emitted gravitational ra- 
diation. Curves are plotted for 10 simulations having different initial black hole 
separations (designated QC-0, etc.) and transition times to perturbative evolu- 
tion. Reprinted from (|9|) and copyright 2002 by the American Physical Society 



(http://link.aps.org/abstract/PRD/v65/el24012). 
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Figure 2: The first gravitational waveform for binary of equal- mass black 
holes evolving through a single plunge orbit, merger, and ringdown, com- 
puted by Pretorius. The waveforms were extracted at four radii from the 
source, and then shifted in time to overlap for comparison. Reprinted with 
permission from (23) and copyright 2005 by the American Physical Society 



( http: / /link.aps.org/abstract /PRL /v95 /e!21 101 ) 
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Figure 3: Puncture trajectories from the merger of an equal-mass nonspinning 
black-hole binary calculated by the UTB group. The apparent horizons of the 
two pre-merger holes, outlined in solid (black) and red (dashed), expand due to 
coordinate effects. Also shown is the first detected common horizon, outlined in 
green (dot-dashed); this designates the point of merger and has a "peanut" shape 
before it settles down. Reprinted with permission from (27) and copyright 2006 by 



the American Physical Society (http://link.aps.org/abstract/PRL/v96/ellll01 ). 
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Figure 4: The universal waveform for equal-mass, nonspinning black holes cal- 
culated by the GSFC group. The figure shows waveforms from four simula- 
tions with increasingly larger initial separations between the black holes. These 
waveforms were shifted in time so that the peak radiation amplitude occurs at 
t = 0. Reprinted from (|34|) and copyright 2006 by the American Physical Society 



( http: / /link.aps.org/abstract /PRD /v73 /el04002 ) 
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Figure 5: Recoil kicks for nonspinning black-hole mergers for mass ratios q € 
{1.0,3.95}, or 7] G {0.161,0.25}; here, rj = q/(l + q) 2 . The open circles are the 
results obtained from numerical relativity simulations by Gonzalez et al. f)61 [) . 
with a functional fit given by the continuous (red) line, and uncertainties of ~ 6% 
indicated by the dashed (blue) lines above and below. For comparison, earlier 
numerical results (discrete symbols with error bars) are also shown. Reprinted 
with permission from (|6ip and copyright 2007 by the American Physical Society 



( http://link.aps.org/abstract /PRL /v98 /e091 101 ) . 
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Figure 6: The longest gravitational waveform for an equal- mass, nonspinning 
black- hole binary merger computed by the Caltech-Cornell group. The left 
panel shows the early stages of the waveform, during the inspiral. The right 
panel displays the merger and ringdown portions of the waveform. Reprinted 
with permission from (46) and copyright 2009 by the American Physical Society 
(|http: / /link.aps.org/abstract /PRD /v79 /e024003|) . 
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Figure 7: Gain in SNR by including the merger segment of the waveform for 
an equal-mass, nonspinning binary. The SNR for sources located at luminos- 
ity distance Dl = 1 Gpc is plotted vs. (redshifted) mass for the Advanced 
LIGO detector. The dashed line shows the SNR calculated using PN tech- 
niques for the early inspiral part of the waveform, — oo < t < — lOOOAf . The 
dotted line shows the SNR using the late inspiral, — 1000M < t < -50M, 
which is the transition region from PN to numerical relativity. The SNR for 
the strong-field merger and ringdown, — 50 M < t < oo, was calculated using 
waveforms from a numerical relativity simulation and is shown using a thin solid 
line. Finally, the SNR from the entire waveform is given as the thick solid 
line.x Here, t = marks the time of maximum gravitational radiation ampli- 
tude. Reprinted from ([73]) and copyright 2007 by the American Physical Society 
flhttp: / /link.aps.org/abstract /PRD /v75 /el24024) . 
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Figure 8: Contours of SNR as a function of redshift z and total binary 
mass M are shown for the LISA detector. These have been calculated us- 
ing a hybrid equal-mass nonspinning waveform: using PN for the early in- 
spiral, matching to a numerical relativity waveform in the late inspiral, and 
continuing with the numerical waveform through the merger and ringdown. 
Reprinted from (|74j) and copyright 2007 by the American Physical Society 
dhttp: / /link.aps.org/abstract /PRD /v75 /e!24024|) . 
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Figure 9: Comparison of merger waveforms for a q = 4 mass ratio non- 
spinning black hole binary calculated using the analytic "effective one-body" 
(EOB) model and using numerical relativity (NR). Reprinted with per- 
mission from (|96p and copyright 2007 by the American Physical Society 
flhttp:/ /link.aps.org/abstract /PRD /v76 /el04049"|> ■ 



